Exact Fidelity and Full Fidelity Statistics in Regular and Chaotic Surroundings 
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For a prepared state exact expressions for the time dependent mean fidelity as well as for the mean 
inverse paricipation ratio are obtained analytically. The distribution function of fidelity in the long 
time limit and of inverse participation ratio are studied numerically and analytically. Surprising 
features like fidelity revival and enhanced non-ergodicity are observed. The role of the coupling 
coefficients and of complexity of background is studied as well. 
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In quantum information it is crucial to know how well 
a prepared state can be isolated, and how the unavoid- 
able mixing with the surrounding behaves. The fidelity 
(often called quantum Loschmidt echo or survival prob- 
ability) describes how the purity of a prepared state de- 
creases due to the interaction with surrounding states 
It serves as a benchmark test of prepared states, 
e.g. a set of qubits in quantum information [2(. Of in- 
terest is to describe the explicit time-dependence of the 
fidelity, and how the decay depends on the coupling to 
the environment, as well as the role played by the com- 
plexity of background states. Time-scales, energy scales 
and explicit shapes of these functions are thus frequently 
studied 0, 0| • Survival probability of a state weakly cou- 
pled to a background has been subject of considerable 
research in mesoscopics and in semiclassics 0, 0, 0, S| • 

F(ermis) G(olden) R(ule) predicts an exponential de- 
cay of the prepared state with decay rate T — 2n\ 2 D, 
where D is the mean level spacing of the background 
and A is the average coupling strength. Deviations from 
this behavior become important when r ~ D [5, 6]. Cor- 
rections to the FGR, which are similar to weak local- 
isation corrections in Quantum Transport can lead to 
non-ergodicity, i. e. the prepared state will never decay 
completely [jj Q. 

The observed saturation of fidelity in the long 
time limit allows us to connect it to the I(nverse) 
P(articipation) R(atio), which is a time independent 
quantity. Usually average behavior is considered. But 
the fidelity in an explicit situation can deviate much from 
the average behavior When constructing quantum in- 
formation and other devices, requiring high fidelity, one 
is interested in a high probability to have states with fi- 
delity superior to some minimal value, above which error 
correction is possible [13]. In this Letter we therefore 
study fluctuations of the fidelity and the full fidelity dis- 
tribution. We find that the latter in the long time limit 
tends to a stable distribution. This distribution is differ- 
ent from the (time-independent) IPR distribution, both 
coincide in the weak coupling limit. In this limit we pro- 
vide an analytical solution. This allows us to study how 
the distribution depends on coupling strength, coupling 



type as well as on dynamics of the background. 

In addition we provide exact solutions for the mean 
fidelity for all times and all coupling strength as well as 
for the mean IPR. We report on nonperturbative features 
such as a fidelity recovery and enhanced non-ergodicity. 

We model the coupling between the prepared state 
with surrounding complexity by the Hamiltonian 

H x = H s + H b + XV = E s \s)(s\ 



N 



+ J2 E v \b v ){b v \ + \J2 (Vu\a){b v \ + h.c.) . (1) 

jy = l !/=l 

where H s represents a special, pure state, that is coupled 
to a background of complex states described by H^, and 
where the coupling, AV, is controlled by the sortless pa- 
rameter A. Without loss of generality we may put the 
unperturbed energy of the special state to zero, E s = 0. 

The Schrodinger equations for the uncoupled Hamilto- 
nians are 

H 8 \s) = and H b \b v ) = E v \b v ) . (2) 
The eigenvalue problem for the coupled Hamiltonian is 



H x \n) = E n \n) 



(3) 



The N + 1 eigenfunctions are expressed in a basis of the 
special state and the background states 



\n) 



,\s) + C nu\K), 



(4) 



(n\s) and c n 



(n\b v 



We model the 



complex surrounding of background states by random 
matrix theory, and describe generic chaotic states with 
an ensemble of Gaussian random matrices that can ei- 
ther show time-reversal invariance (GOE, (3 — 1) or not 
(GUE, (3 = 2). Regularity of the background states is 
modeled by assuming Poisson statistics for In all 
cases the spectrum is unfolded so the mean level spacing 
equals one, D = 1, at least in a surrounding of the spe- 
cial state. This implies that the energy scale, including 
the coupling strength, is always expressed in units of the 
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mean level spacing, while the time scale is expressed in 
units of the Heisenberg time, th — h/D. In the following 

n = i. 

Matrix elements of the operator V between the spe- 
cial state and the complex surrounding, V v = (s|V|6„), 
are taken as Gaussian distributed random numbers with 
zero mean and variance one. As we will see, it is impor- 
tant to distinguish between real coupling (V v € K) and 
complex coupling (V„ € C). The size N of the Hilbcrt 
space describing the complex surrounding is in principle 
infinite in RMT. In the numerical simulations it is taken 
sufficiently large to achieve convergence. 

We first focus on fidelity decay. The fidelity ampli- 
tude f\(t) = (&o(t)\^!\{t)) is the overlap, between the 
pure state, time developed under the influence of the 
pure Hamiltonian, |*o(*)> = exp (-i(H s + H b )t) 1*0(0)), 
and under the influence of the perturbed Hamiltonian, 
\ty\{t)) = exp(—iH\t) |*o(0)). It is a measure of how 
the pure initial state, |\&(0)), gets disturbed or mixed 
due to the (unavoidable) coupling to the complex sur- 
rounding a time t later. We are interested in the decay 
of the special state \s), i. e. 1^(0)) = \s). We expand it 
in a basis of eigenstates to H\. The fidelity amplitude 
becomes 



= ^2\cns\ 2 - exp (iE n t) 



(5) 



It is seen that f\ (t) is the Fourier transform of the local 
density of states (LDOS), 



p{E) =Y,Ks\ 2 6(E-E n ). 



(6) 



The smooth part of the LDOS follows a Breit-Wigner dis- 
tribution, with the width given by T — 2ttX 2 , as obtained 
from FGR under very general assumptions. Therefore 
the mean fidelity amplitude of the special initial state | s) 
will unavoidably decay exponentially 



A(t) = exp(-rt/2) , 



(7) 



where the bar denotes average over background and over 
coupling matrix elements. 

The fidelity (survival probability), of the special state 
is defined as, F\(t) — \f\(t) | 2 . In a Drude-type approxi- 
mation for the mean fidelity 



/a =exp(-R) 



(8) 



FGR is recovered. This result is also obtained by second 
order perturbation theory. We write fidelity as the sum 



F\(t) = IPR + Ff\ uc (t) , 



(9) 



of a constant term and a term which is fluctuating on a 
timescale comparable to Heisenberg time. The constant 
term 



IPR = ^|c ns | 4 = £ JdEp 2 (E), 



(10) 



is the inverse participation ratio of the special state in 
the basis of the eigenvectors of the full Hamiltonian. The 
fluctuating term 

F auc (t) = 2j2Ks\ 2 \c ms \ 2 cos((E n - E m )t) (11) 



vanishes, if we average fidelity over a time window, 
large compared with Heisenberg time. Using a Breit- 
Wigner distribution for p(E) and the Drude approxima- 
tion Eq. ©, IPR7 = D/(ttT) is found @], which can 
obviously not hold for small T. 



For complex coupling we were able to calculate F\(t) 
exactly for a regular and for a GOE/GUE background. 
For a regular background the result ll| is given by 



Fx(t) =1- 



A f dx _i 



J V% 



_ e 4(1-*) 
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2t t 2 V> 



e x cosh f 



- 1 



sinh 



s/T 



(12) 



For a GUE/GOE background the corresponding more 
complicated expressions can be found in [111 ]. For small 
times this function decays exponentially according to the 
FGR law. Surprisingly, after some characteristic time 
fidelity reaches a minimum and increases afterwards to a 
A dependent saturation value F\(oo) = IPRa. We find 
for a regular background 



IPRa = 1 



/tF 3 A d /ttA 
2 V~2~ 



(13) 



with D(w) = exp(w 2 )erfc(a;). IPR A is a monotonously 
decreasing function. For small coupling the saturation 
value behaves as IPRa — 1 — 7r 3 / 2 A. For large values of A 
it decays algebraically as ~ 2/(7r 2 A 2 ) which is four times 
the value, obtained by the Drude approximation, Eq. (jSJ. 
This is a striking enhancement of non-ergodicity. The 
corresponding expresssion for a GUE background is 



Fa(oo) = 1-vr^ 



^ 3 X 



1 - 2tt 2 A 2 ) D(ttA) (14) 



This function decays algebraically as ~ l/(7r 2 A 2 ) for large 
A which is twice the value, predicted by the Drude ap- 
proximation. 

In Fig [1] fidelity F\(t) is plotted for A = 0.1 for different 
complexity of the background. For complex coupling the 
curves are obtained from Eq. (|12p respectively from the 



corresponding expressions taken from Ref. 11(. For real 
coupling the curves are obtained by Monte Carlo simu- 
lations. It is seen that in all three cases fidelity reaches 
a minimum and saturates afterwards at a finite value. 
Fidelity decay is stronger for a chaotic background than 
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FIG. 1: (Color online): Evolution in time of fidelity for cou- 
pling strength A = 0.1. The full lines describe top down 
fidelity decay for real coupling to a Poissonian (green), to a 
GOE (red) and to a GUE (black) background. The dashed 
lines describe top down fidelity decay for complex coupling to 
a Poissonian (green), to a GOE (red) and to a GUE (blue) 
background. 



for a regular background. This is in accordance with the 
original perturbative arguments by Peres 

Quite remarkably the decay of fidelity is much less sen- 
sitive to the complexity of the background than to the 
structure of the coupling. There is practically no dif- 
ference between a time reversal invariant chaotic back- 
ground and a background with broken time-reversal in- 
variance. Nevertheless the difference between a real cou- 
pling of the special state to the background and a cou- 
pling which breaks time reversal symmetry is sizeable. 
For one reason, because the return probability from the 
background into the special state is suppressed by a cou- 
pling, which breaks time reversal invariance. 

We now turn to full fidelity statistics. We introduce 
the full distribution function 



P F (c,t) = 5(c-F x (t)) 



(15) 



The distribution of the fidelity Pf (c, t) is calculated nu- 
merically and plotted for different times in Fig. [2] A 
saturation of the distribution function is found for times 
larger than a (A-dependent) saturation time. In Fig. [2] 
this saturation time is about 20 times Heisenberg time. 
The saturated distribution is shown in Fig. [3] for GOE 
statistics and for coupling strength A = 0.05. It can be 
compared to the IPR-distribution 



Pipr(c) = S (c-J2\ 



(16) 



We see that the two distributions are similar but dif- 
ferent. The reason lies in the fluctuating term Ffl uc 
(see Eq. ([9])). Although the ensemble average of Ffl uc 
vanishes, its variance does not. Ultimately Ffi uc (t) con- 
tributes substantially to the full fidelity distribution. As 



FIG. 2: (Color online): Fidelity distributions at times t = 1 
(pink), t = 2.5 (blue), t = 5 (red), t = 10 (green) and t = 20 
(black) for coupling strength A = 0.1. For times t > 20 the 
distribution is stable. 



a result the stable fidelity distribution and IPR distribu- 
tion are different. 

We are interested mainly in cases of high fidelity, that 
is, when the coupling strength, A, is small. In this limit 
Pf(c, oo) should be better and better approximated by 
Pipr(c). In this limit an analytic result for Pipr(c) can 
be obtained. By solving the Schrodinger equations ([2]) 
and ([3]) an expression for the component of the special 
state in the eigenstates \n) of H\ is obtained as 



N 



= h^ 2 E 



IK 



J (E n - E^f 



(17) 



This exact expression for the components contains eigen- 
value solutions, E n , as well as input matrix elements V v 
and energies E v . For small A we may approximate the 
exact eigenvalues E n by their unperturbed value E v t n )- 
If we denote by |0) the eigenstatc of H\ which has evolved 
from the special state | s) , we see that in this approxima- 
tion all |c ras | but |co s | vanish. This means that for small 
couplings the IPR will be dominated by only one term 

I COs 



I 4 . We define the distribution 



P (c) = 8{c- |c 0s | 4 ) . 



(18) 



Then for small A, Pf(c, oo) ~ Pipr and likewise Pipr 
~ Pq. We could perform the ensemble average of Pq 
exactly. For detail of the calculation see 12j. For a 
regular surrounding (Poisson statistics) we find: 



1 



2\a f 



4^(1-V^) 3/2 



(19) 



where a\ — for real coupling and a2 = yft/i for 

complex coupling. For chaotic surrounding we find for 
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FIG. 3: (Color online): Comparison of the distributions Po 
obtained from Eq. (|20|l (dashed red line) with Pipr (full green 
line) and Pf(c, oo) (full black line) obtained from Monte- 
Carlo simulations for a coupling constant A = 0.05. 



time-reversal symmetry (GOE): 



Po(c) 



n 3 \ 6 c 



-x 



4^3 V 2 (!-^) 5 



1 [KoiX!) + Ki(X!)) 
(20) 

where K n are modified Bessel functions of second kind. 
For a chaotic background without time reversal symme- 
try (GUE) we get: 



Po(c) 




2ttA 2 



-X-2 



1 + 2A 2 ), (21) 



with Xr- 



V P ~ ^i{i-y/c) ' ^ n El we compare the ana- 
lytical expressions for the distribution Po to numerical 
simulations of Pipr, and of Pf(c, oo) for A = 0.05 in the 
case of a GOE background. We observe good agreement 
of the analytically calculated Po (c) with the distribution 
Pipr(c) obtained by simulations. Both curves are in- 
distinguishable at least in the range c > 0.6. There is a 
small but notable difference to the distribution Pf(c, oo), 
which vanishes if we go to smaller values of A. 

In conclusion, we studied fidelity decay for a special 
state coupled to a regular or chaotic environment. We 
found a saturation in the long time limit and a revival. 
In [l]| a fidelity freeze was predicted. For a purely 
off-diagonal perturbation, after an initial decay fidelity 
freezes on a plateau for some time and decays afterwards 
to zero. In the model considered here, the perturbation V 
is purely off-diagonal as well. The saturation, we found 
here might thus be considered an extreme case of fidelity 
freeze. The fidelity revival found here is genuinely dif- 
ferent to the one reported earlier 14|. There, a satisfac- 
tory explanation was given by the spectral rigidity of the 
GUE/GOE0, HH. The fact that the revival occurs for 
a regular background as well encumbers such an expla- 
nation in the present case. 



The saturation of fidelity is a direct consequence of the 
fact that the fidelity distribution relaxes in the long time 
limit into a stable distribution. In the small coupling 
limit it becomes the distribution of the IPR. In this limit 
we found an analytic expression. Both distributions have 
a rich structure and are highly sensitive to even small 
changes in the coupling strength. Their relation to max- 
imum strength distribution, introduced recently 3, 17 1 



will be discussed elsewhere 121. Their calculation for 



arbitrary coupling strength is a challenge for the future. 

Our results yield an important benchmark for the de- 
cay of a prepared quantum state. It might be probed 
numerically and experimentally in chaotic quantum sys- 
tems [l8, lj| 2(| or on quantum information devices. One 
instance for a possible numerical experiment is the decay 
rate of a regular state in a mushroom billiard due to 
dynamical tunnelling into the chaotic part of the phase 
space. 
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